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FOREWORD 


This report presents the results of work done by Northrop Services, Inc., 
Huntsville, Alabama, for the National Aeronautics and Space Administration, 
George C. Marshall Space Flight Center, under Contract NAS8-21810. The work 
was performed for the Science Aid Engineering Directorate in response tc> 
Appendix A, Schedule Order Humber A02Z(A-13) 

Dr. George H. Fichtl was the Technical Coordinator for thib task. The 
author is grateful to Dr. Fichtl for his guidance and fruitful discusnions 
during this work. 
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Section I 
INTRODUCTION 


Future studies of the ascent control of the Space Shuttle are being 
planned. Some of these planned studies will be carried out analytically using 
a model simulation procedure. In this technique the Space Shuttle is followed 
analytically as it ascends through the atmosphere. The model interacts with 
the simulated random winds, and the ability of the ascent control to maintain 
the desired trajectory is studied. 

The simulated winds must be generated so as to have the appropriate 
statistical behavior. This report discusses the method of generating a 
random wind signal using a digital computer that will have the same statistics 
as the winds encountered in the Space Shuttle ascent. 

The procedure for generating the random wind is to develop a control 
system which inputs discrete white Gaussian noise and outputs a random signal 
that has the statistical behavior of the wind. The control system is written 
in terms of state equations which are then digitized for computer calculations. 


Section II 

SIGNAL ANALYSIS 


2.1 EMPIRICAL AUTOCORRELATION 

An empirical wind autocorrelation was obtained from detailed Jimsphere 
measurements (Reference 1) made at the Kennedy Space Center and shown on 
Figure 2-1. The random components of the longitudinal winds, u and v, were 
obtained as a function of altitude z. The turbulent wind component was 
normalized to yield 

»«*> ■ 5 $ 

where t ■ z/L(z). The term L(z) is a length scale which is chosen so that the 
dimensionless y process is homogeneous; that is, the second order statistics of 
y are independent of z. The o(z) is the standard deviation of v(z). 

The resulting empirical autocorrelation is given by 

R y ( t l» t 2 ) - <y(tj) y(t 2 )> - <y(tj) y(t x +T)> - R y (r) 

where <•> represents an ensemble average and t is the lag. The same empirics : 
autocorrelation was found to apply to both the u and the v components of thu 
horizontal wind when they were appropriately normalized (Reference 1). 


The empirical autocorrelation can be represented in a functional form 
which can be Fourier transformed to give the power spectrum for the dimension- 
less wind. By factoring the power spectrum, a control system function can 
be obtained. A control system is then defined which Inputs white noise and 
outputs a random signal which has the same autocorrelation as the functional 
autocorrelation just discussed. The control system can then be written in 
terms of stats equations which are put in a discrete form for use on a digital 
computer. The autocorrelation of the digitized output signal is in good 
agreement with the desired autocorrelation. 
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Figure 2*1. AUTOCORRELATION OF WIND SIGNAL 





2.2 FITTING EMPIRICAL AUTOCORRELATION 

An empirical autocorrelation was obtained from wind data (Reference 1) 
and la plotted on Figure 2-1. The empirical autocorrelation waa approximated 
In functional form by 

K y (t) - <y(t) yCc + x)> 

- exp (-d|t|) { coa B(t) - g ain b|t|} (2-2a) 

where the notation <*> ref era to eneemble averagea and B and 0 are empirically 
determined coefficienta. 

* ■ * <2-2b) 

D - 0.539 

and t ie the autocorrelation lag. A coaparlaon of the empirical and functional 
form of the autocorrelation In Figure 2-1 ahows good agreement between the two. 


By taking the Fourier tranaform of the autocorrelation the power apect-um 
♦y la obtained; and, aince By la an evan function, the power ape ct rum can be 
written as 




L\ 


(t) coa WTdT - 


4&j£ 

[ D 2 + (B-ui) 2 ] [d 2 + (B + w) 2 


(2-3) 


2.3 SYSTEM FUNCTION 

In general, the output power epectrua of a ay a tea can be written aa 

♦ y ■ Fj (2-4) 

where la the input power apactrum and H la the ayatem function. The H* ia 
the complex conjugate of the ayatem function. 
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If the input le white Gaueeien noise, then ♦j « 1.0; and $ y equals HH* 
where H is given by 

„ (2 Jpj S 


(8 + D - IB) (S + D + iB) 3 

where S • iw. This result is shown by Figure 2“2 . The systes equation can 
now be written as 

Y(S) - H(S) 1(8) 

where Y(S) is the diaensionless wind having the autocorrelation given by 
Equation 2-2 and 1(8) is the Gaussian white noise Input (Figure 2-3). 


i i 1 1 mi 


Figure 2-2. MAGNITUDE OF SYSTEM FUNCTION 
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2 A STATE SPACE SYSTEM 

As shorn in Reference 2, m csn, dsflns ths system In Corns of ststo 
vsrisbls X t givon by the following equstion where the Einstein suanetion 
convention is inplled by repented indices: 



\j + d l 1 


(2-7) 


The eye ten output is given by 

T - e A X ± (2-8) 

Following Dorf (Reference 2) the flow grsph stste nodal is produced as shown 
on Figure 2-4, where 

b. - 2 /IT; • - 2D *nd • - D 2 + B 2 (2-9) 

. 1*0 
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Section III 

DISCRETE TIME SYSTEM 


r 


3.1 DISCRETE STATE SPACE SYSTEM 

For use on digital coaputers the state equations aust be converted to a 
discrete tiae systea. One procedure for achieving this is given in Reference 
3. In this procedure the input signal I ia passed through a zero order 
holding device which aaaples the signal at unit intervals of tiae and holds 
the signal value constant between saaples (Figure 3-1). The procedure for 
converting to a discrete tiae eystea is discussed next. 



Figure 3-1. CONTROL SYSTEM WITH SAMPLER AND HOLDING DEVICE 


Equation 2-7 can be integrated as shown in Reference 4 to give 


*i(t) • ♦jj (* * l # ) Xj (* 0 ) + " T ) d|I(t) 


dr 


(3-1) 


Miere is known as the fundamental aatrix. Since I is considered constant 
over the interval T, Equation 3-1 can be evaluated at tins t ■ (K ♦ 1)T 
over the increment T and obtain 


* 4 (E ♦ l) - (T) Xj (X) ♦ A t (T) I(X) 


(3-2) 
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\ 




where 



<t) dj dt 


( 3 - 3 ) 


Following the usual procedures, the fundaaentel natrix can be evaluated using 
Laplace transforas where L represents ths Laplace transfora operation. 



where 


. JO if i ¥ J 
iJ [lifi-J 


( 3 - 4 ) 


This can be solved to give 
-DT 




t 


e __ | cos BT — sin VT 


■] ‘ 


e _DT sin IT 


(3-5) 


- y- e" OT sin IT 


-DT 


cos BT - | 


sin BTl 


Taking the Halt for saall T 


E 


♦u«>l - 


•si 1-2DT 
o 


M 


• o<T«l 


( 3 - 6 ) 


Then, froa Equation 3-3, obtain 


h-i: C-J " ■ 0 ■ [°J ! 

This then gives the result for ths discrets case as 


o<T«l 


(3-7) 


X t (K + 1) - A 4J Xj(K) ♦ D 1 I(K) 

Where A^j and are defined by Equations 3-6 end 3-7. 

* 


(3-8) 
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This same relationship can be obtained by using a forward finite 
differencing technique (Reference 2). In this method Equation 2-7 can be 
written as 

X <K + 1) - X (K) 

— * a tJ Xj(K) + d t I(K) (3-9) 

This expression can be rewritten to give the same result as Equation 3-8. 


\ (K + 1) - (Ta iJ + & u ) Xj(K) + Td ± I(K) 
- k u Xj(K) + D t I(K) 


(3-10) 


3.2 EFFECT OF DIGITIZING ON AUTOCORRELATION 

As shown on Figure 3-1, a zero order holding device has been added 

(Reference 4). Therefore, the input into the continuous system will be put 
in a discrete form consisting of a stair junction of Gaussian heights and of 
width T (Figure 3-1). As discussed in Reference 5 a correction must be made 
for this effect. This can be done by finding the spectrum of the discrete 
input. The autocorrelation of the discrete input is given by 


(t) - <1* (t) I'(t + t)> 
- a 2 Pr[A] 


(3-11) 


2 

The a is the variance of the Gaussian noise input, and Pr[A] is the probability 

that points t and t + t of the input in the discrete form both occur between the 

times KT and (K + 1)T. Since the Gaussian noise input has a variance of 1.0 
2 

then o » 1. The probability of \ occurring (Reference 6) can be seen to be 


Pr(A) 


f 

i 


T- 111 

T 


|t|<t 

M>t 


Then the input autocorrelation le given by 


(3-12) 


Ri<t) 



(3-13) 
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This can he Fourier transformed to give the Input power spectrum 


♦1 


aln (u»T/2) l 
Tw 

2 J 


' T ; o<T<<1 


(3-14) 


Thla same result was found In Reference 5. 


Thus , the output spectrum of the discrete input, y ' , is now given by 


♦ y i - THH* - T 4 y 


(3-15) 


Therefore, the digitised autocorrelation is 

R y ,(x) “ TR y (t) (3-16) 

Thus, the output signal must be normalised to obtain the correct autocorrelation 
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Section IV 

THEORETICAL DISCRETE AUTOCORRELATIONS 


The theoretical autocorrelation can be calculated for the discrete 
equation following the procedure in Reference 5 by writing Equation 3-8 as 


X t (K + 1) - A tJ X^K) + Dj I (K) 


X 1 (K + 2) - A aj Xj(K + 1) + D i I(K + 1) 


A iic V K > + A iK °K 1(k) + D i I(K + 1} 


(4-1) 


where 


A^« A A 
IK iJ JK 


(4-2) 


Similarly 


(K + 3) 


A iK >X K< K) + A iK >D K 1(K) + A iK D X I(K + X) 


+ D i I(K + 2) 


(4-3) 


and, in general. 


K+n-l 


(K ♦ n) ■ + I A 1J (K+ "- l - r> D J I(r) 


(4-4) 


where 


*U >D J ■ D i 


(4-5) 


Letting K + n - i, and, as R+-«, X (-•) - 0 


l-l (i-l-r) 

X t (t) • l K u Djl(r) 


(4-6) 
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This gives a nonrecursive form for obtaining (l). If m - t-l-r then the 
equation becomes 


X i U) - l A i j" ) D I(l-m-l) 
m-0 


(4-7) 


Since <I> ■ 0 then <Xj> ■ 0 


The autocorrelation can be found as follows: 
using Equation 4-7 


<1(0 X.(l+X)> - R' (X) - l A. j D T <I(Jt) I(i + X - m - 1)> 
1 1 l m-0 1J J 


(4-8) 


Since I is Gaussian white noise 


<I(l)I(4+X-m-l)> - 6(X-m-l) ; where 6 


So that Equation 4-8 becomes 


0 when X-m-1 ^ 0 


1 when X-m-1 - 0 


(4-9) 




f ° r 

0 for X-0 


Similarly, 


^(X) - <X i (i) X K (i+X)> - [ D J <I(i-m-l) X R (l+X)> (4-11) 


J 0 ‘J’ 1 R lX K < X+ * tl > 


*i.x< x) ’ A A u’ > D J A xt + * > D » ' l f <"> *<•«> 


(4-12) 


The result can be seen to be a convolution suaeutlon which is the discrete 
time counterpart to the convolution Integral. 
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For Che autocorrelation of the system output y 

M1) . o»> >(»»). . "'y .. 


,2 


.2 


(A-13) 


Then, by combining Equations 4-12 and 4-10, obtain 

.2 m 2 “ # (m) A (mfX) 
m-0 


u* -.2 r 

D 1 1 4 a 22 a 22 


I A, 


(m) .(ort-A) 
22 22 


R (A) - 
7 


,2 


l 1.4? ) 2 


m-0 


(4-14) 


where the normalizing factor is given by 


°' 2 * <y 2 (Jt)> - b* Rj.^ -b^T 2 l (A^) 

y 1 *2*2 1 m-0 


(4-15) 


The analytical result for the discrete autocorrelation R^(A) is shown in 

Figure 2-1 for T • 0.125 and T - 0.1 and is in good agreement with the desired 

2 

autocorrelation. The values found for o' were 0.145 for T of 0.125. 

2 * 

For T - 0.1, 0 ^ * 0.111. This indicates that the system is stable since 

the variance is finite. The correction factor found earlier in Equation 3-16 
2 

was - T, which in the present case gives close agreement with the more 
detailed discrete results. 
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Section V 

STABILITY ANALYSIS 


In Reference 4 it is shown that a stationary linear system subjected to 
a bounded input is stable if and only if all the zeros of the characteristic 
polynomial, “ ^jl within the circle |x| * 1 in the complex X plain. 
This results in 


XI 


6 J 



(5-1) 


where | * | denotes a determinant . The eq • 5-1 can be expanded to 
X 2 + X(2DT - 2) + 1 - 2DT + a Q T 2 - 0 


(5-2) 


This can be solved to give 

X - (1 - DT) + iTB (5-3) 

which has the magnitude 

| X | - Tl - 2DT + T 2 (D 2 + B 2 )j 1/2 (5-4) 

Since (D 2 + B 2 )T < -2D for small T then |x| <1 ; and therefore, the system 
is stable. 
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Section VI 


COMPUTER SIGNAL OUTPUT 


Discrete Gaussian white noise can be generated on the computer using 
readily available programs. Inputting this into the recursive equation for 
(Equation 3-8) can result in a set of discrete values of y(K); K ** 1,2,3... 
This result can be normalized either by o^ (Equation 4-14) or by given by 



The resulting autocorrelation was calculated from 


R (A) 
y 


h I y(K) y(K+A) 
k-l 


For a time Increment of T ■ 0.125 and 1,000 samples the result in Figure 2-1 

was obtained. This result is in good agreement with the desired autocorrelation. 

2 

The value obtained for for T ■ 0.125 waa 0.142 which is in good agreement 
with the result obtained theoretically which is given in Section IV as 0.145. 
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Section VII 
CONCLUSIONS 


A technique is developed for the s Wation of random wind signals having 
appropriate autocorrelation which can be readily generated on a digital 
computer. These results are to be used in generating wind data tapes for 
Space Shuttle launch simulations. These results can be linearly interpolated 
to give intermediate values between the generated results. 
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